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ABSTRACT 

We present an analytic formalism that describes the evolution of the stellar, gas, and 
metal content of galaxies. It is based on the idea, inspired by hydrodynamic simula- 
tions, that galaxies live in a slowly-evolving equilibrium between inflow, outflow, and 
star formation. We argue that this formalism broadly captures the behavior of galaxy 
properties evolving in simulations. The resulting equilibrium equations for the star 
formation rate, gas fraction, and metallicity depend on three key free parameters that 
represent ejective feedback, preventive feedback, and re-accretion of ejected material. 
We schematically describe how these parameters are constrained by models and obser- 
vations. Galaxies perturbed off the equilibrium relations owing to inflow stochasticity 
tend to be driven back towards equilibrium, such that deviations in star formation rate 
at a given mass are correlated with gas fraction and anti-correlated with metallicity. 
After an early gas accumulation epoch, quiescently star-forming galaxies are expected 
to be in equilibrium over most of cosmic time. The equilibrium model provides a sim- 
ple intuitive framework for understanding the cosmic evolution of galaxy properties, 
and centrally features the cycle of baryons between galaxies and surrounding gas as 
the driver of galaxy growth. 



1 INTRODUCTION 



Galaxy formation involves a wide range of diverse physi- 
cal processes operating on stellar to cosmological scales, in- 
cluding the hierarchical growth of structure, star formation, 
black hole accretion, and a plethora of poorly-understood 
feedback processes that strongly modulate galaxy growth. 
Given this complexity, it is surprising that galaxies display 
simple and tight scaling relations between many of their 
key constituents. These include well-established relations be- 
tween the bulge velocity dispe r sion and central bla ck hole 
mass (e.g. iGultekin et al.ll2009l: iGraha m et aLl lgOllI ). circu- 
lar velocity and luminosity ( TuUy fc Fisher 19 771), star for- 
mation rate and stellax mass (e.g. lDave 2008; Gonzal ez et al.l 
[2010, and references therein), and m etallicity and stellar 
mass (e.g. iTremonti et all [2004 : lErb et all 120061 ). Each has 
low scatter and evolves roughly independently of mass. The 
simplicity of these relations hints at an underlying unifor- 
mity in galaxy evolution that is not immediately evident 
from the complexity of current hierarchical galaxy forma- 
tion models. 

In the longstanding canonical scenario for galaxy forma- 
tion, galaxies form as angular momentum-conserving disks 
cooling from hot gas bound within dark matter halos, and 
these disks su bsequently merge to fo r m larger and earlier- 
type galaxies llRees fc Ostriked Il977l: I White fc Reed [l97i : 



IWhite fc Frenkill99ll : IMo. Mao, fc White! Il998l ). This sce- 



nario is well-situated within currently favored hierarchical 
cosmologies, and analytic models based on it have been quite 
successful at repr oducing man y observed galaxy proper- 
ties (see review bv lBensonll2010l ). However, the present gen- 
eration of such models (often called "semi-analytic" mod- 
els) are typically enormously complex, with a host of free 
parameters describing various interrelated physical phenom- 
ena. Numerical simulations that explicitly track gas dynam- 
ical processes enable a more ab initio calculation, but still 
require many "sub-grid" parameters for key physical pro- 
cesses and are in practice limited by dynamic range and nu- 
merical uncertainties. In either case, the complexity of such 
models makes it difficult to extract simple physical intuition 
for what drives the evolution of basic galaxy properties. 

In this paper, we present an analytic framework for un- 
derstanding the evolution of the stellar, gas, and metal con- 
tent of galaxies. This framework is based on intuition gained 
from hydrodynamic simulations of galaxy formation. In such 
models, galaxies are fed primarily by cold (~ 10 ^ K) streams 
conn e cting to filamentar y large-scale structure (|Keres et al.l 
I2OO5I : iDekel et al.l |2009|). ou tflows are strong and ubiq- 
uitou s ( Springel fc Hernguistl l2003bl : lOppenheimer fc Pavel 
|2008| ). and outflowing material commonly return s to galajc- 
ies ( "wind recycling" ; lOppenheimer et al.ll201Gl ). Hence in 
this framework, galaxy evolution is governed by the cycle 
of baryons exchanging matter and energy between galaxies 
and surrounding intergalactic gas. 
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Our framework is an attempt to distill the insights 
gained from such hydrodynamic simulations into an analytic 
formalism that both describes the results of simulations and 
provides intuition into the key physical d rivers. It is based 
prima rily on the formalism presented in iFinlator fc Davj 
lj2008h . with key ext ensions, and shares featu res with vari- 
ous recent works (e. g. Rasera Sz T cvssicr 2006; Bouche et al.l 



I2OIOI : iDutton et al]|2010l : lKrumholz fc Dekelll2011 V indicat 
ing a groundswell towards this "baryon cycling'lj view of 
galaxy formation. We demonstrate that the star formation 
rate, gas content, and metallicity of galaxies can be de- 
scribed by simple equations that depend on three param- 
eters that are directly related to inflows, outflows, and wind 
recycling. These parameters are poorly constrained in both 
value and functional form, and hence the number of free 
parameters in this model may be much greater than three, 
pending observations that can better constrain them. To- 
gether, these parameters quantify the impact of baryon cy- 
cling on galaxy growth. We give examples of how these equa- 
tions lead to straightforward intuitive explanations, often 
differing from traditional ones, for the results seen in recent 
observations and models of galaxy evolution. 

This paper begins in 311 by describing the basis for our 
analytic framework, namely the equilibrium condition, and 
discusses the physical constraints on its various terms that 
ultimately govern stellar growth. In ij3] we present an ex- 
pression for gas fractions and explore some implications. 
31] discusses what governs galax;y metallicities, and relates 
this to wind recycling. ^S] gives some brief examples of how 
these equilibrium relations yield straightforward intuition 
into what governs basic galaxy properties. SJS] discusses a 
preliminary implementation of an equilibrium model, and 
explores some parameter variations. i}7] considers what hap- 
pens when galaxies depart from equilibrium owing both to 
stochastic fluctuations and more permanent departures. 31] 
discusses when galaxies first attain equilibrium in the early 
universe. Finally, we summarize and discuss broader impli- 
cations of our framework in 33 



2 THE EQUILIBRIUM CONDITION 

Star-forming galaxies in hydrodynamic simulations are usu- 
ally seen t o lie near the equilibriu m co ndition (see e.g. Fig- 
ure 13 of iFinlator fc Da^ |2008| . and iDutton et alTbOlOl : 
iBouche et al.ll2010l ): 



(1) 



where the terms are the mass inflow rate, mass outflow rate, 
and star formation rate (SFR) , respectively. Inflow and out- 
flow refer to gas motion in and out of the galaxy's star- 
forming region, i.e. the interstellar medium (ISM). Star- 
forming galaxies fluctuate around this relation but are gen- 
erally driven back to it on short timescales, as we discuss in 
37] this "self-regulating" behavior is why we dub this model 
the equilibrium model. 

The equilibrium condition is close to an expression for 
mass conservation, except that it importantly does not con- 
tain a term describing a gas reservoir. A key ansatz of this 



^ A term coined in the Astro2010 Decadal Survey Report: New 
Worlds, New Horizons 



formalism is that the rate of change in the gas reservoir is 
small compared to the other terms in Equation [T] The mo - 
tivation for this ansatz is that in IFinlator fc Pavel (|2008l ). 
we found this to be explicitly true in hydrodynamic sim- 
ulations of galaxy formation. We note that this scenario 
has als o been referred to as a "reservoir" or "bat h tub" 
model l|Bouche et al.ll2010l : iKrumholz fc Dekelll201lh . 

Defining the mass loading factor rj = Afout /M« , we can 
rewrite the equilibrium condition as 



SFR = Min/(l-|-r?). 



(2) 



Hence in this scenario, a galaxy's star formation history over 
cosmic timescales is determined by the evolution of A/in and 

Let us consider infiow first. Afin can be broadly sepa- 
rated into three terms: 



A'4, 



Baryonic inflow into galaxy's halo, which is 



primarily set by the a ssumed cosmology . Here we employ 
the form forwarded bv lDekel etlo] (120091 '): 



= 0.47/5 



0.15/1 + 2 



A^halo V 1012 Af© 

iFakhouri. Ma. fc Bovlan-KilchinI (|2010l ') presented 
a different parameterization of 25.3M^;^^{1 + 

1.65z)x/n^(l -f z)3 -I- Mq/jt, while 

iFaucher-Giguere. Keres. fc Mai (|201lD found 33.6M^^{1 + 
0.91z)^r2m(l -I- z)3 -I- JIa Mq/jt; these yield similar results 
over most of cosmic time. 

• A/prov = The amount of the gas entering the halo that 
is prevented from reaching the ISM (to be subtracted from 
the halo infall). This is material that ends up in the gaseous 
halo of the galaxy, and can also be regarded as the rate 
of growth of halo gas. We characterize this by defining a 
preventive feedback parameter 



Afprcv/A/g: 



(4) 



Gas infalling that has previously been ejected 



in outflows, along with gas returned to the ISM via stellar 
evolution. This provides an extra component in addition to 
the baryons associated with the gravitational infall of dark 
matter (i.e. A/grav). 



Expressing Afin in these terms, we obtain 



Min = Afg: 



Afp 



+ Mr, 



CAfgrav + A/r< 



(5) 



Halo infall (Afgrav) is driven b y gravity, hence is mostly 
independent of feedback process es ( van de Voort et al ]|201ll : 
iFaucher-Giguere. Keres. fc Mai |2011^ and is determined 
primarily by cosmology and the halo merger rate (e.g. 
iNeistein. van den Bosch, fc Dekel 1 12OO6I: INeistein fc Dekel I 



bood lMcBride. Fakhouri. fc Mall2009l : iGenel et al.ll2O10D . Iii 

contrast, Afprev and A/recyc are direct consequences of feed- 
back processes. We will discuss Afrecyc further in 33] when 
we will relate it to the metallicity of the infalling gas. This 
leaves the preventive feedback parameter ij, which we now 
consider. 

There are a number of sources of preventive feedback, 
each with its own dependence on halo mass and redshift, the 
products of which comprise the total 

• Cphoto represents suppression of inflow owing to photo- 
ionisation heating. This operates at low masses, and ap- 
proaches zero below a photo-suppression mass that increases 
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Figure 1. Schematic illustration of the preventive feedback pa- 
rameter C(-'^^halo)- Cphoto (magenta) is parameterized here as (1 + 
UM /M-t)'^)-'^-^ (following eq. 1 of (okamoto. Gao. fc Theunj 
1 20081) . and Cquench (red) is analogously parameterized. Cgrav 
(green) is taken from Equation (6] Cwinds (blue) is an arbitrary 
function depicting an increasing effect to smaller masses from 
wind heating. Top panel: We choose M^^ = IO^Mq and Mq = 
10^^ Mq, perhaps appropriate for z = 2. Bottom panel: We choose 
Mj = 3 X 10^ M0, Afq = IO^^-^Mq, and a stronger Cwindai ^ore 
appropriate for 2 = 0. These values and their parameterizations 
are fairly arbitrary, because the terms are poorly constrained; 
this plot is purely intended to illustrate general trends. The total 
C = Cphoto Cquench Cgrav Cwinds is shown as the black line, showing 
that vigorous star formation can only proceed within an interme- 
diate range of halo masses where approaches unity. 



from a halo mass of ~ 10 Mq during reionisat i on to 
few X 10^ Mq at the pr esent epoch (|Gnedinl [2OO0I : 



lOkamoto. Gao. fc Thetmil2008l ). 

• Cquench IS assoclated with whatever physical process (es) 
quench star formation in massive halos and prevents cool- 
ing flows, probab ly related to feedback from supermassive 
black holes (e.g. ISomerville et al.]|2008l ). It drops to zero 
abov e the quenching mas s Mq ~ IO^^Mq (e.g. lCroton et al.l 



I2OO6I: iGabor et al 



■ 201 ll ). which may be higher at high- 
z (|Dekel et all |2009| ). Th is may further dep end on the 
merger history of galaxies (|Hopkins et al.|[2008l ). 

• Cgrav reflects suppression of inflow by ambient gas heat- 
ing owing to gravitational structure formation via the for- 



mation of virial shocks. Using hydrodynamic s i mulat ions 
with no outflows, iFaucher-Giguere. Keres. fc Mai l|201ll ) de- 
termined 



Cg: 



0.47 



l + z 



7 Mh 



V1012M, 



(6) 



We note that those simulations did not include metal- 
line cooling, which may be an important effe c t for 
heating gas in virial shocks (|Dekel fc BirnboimI I2OO6I : 
lOcvirk. Pichon. fc T evssier' 2008); nonetheless, we show in 
Figure [2] that our simulations including metal-line cooling 
yield similar results. 

• Cwinds is associated with additional heating of sur- 
rounding gas provided by energetic input from winds. This 
tends to affect lower-mass systems more, but is highly 
dependent on the physics of how outflows interact with 
surrou nding gas, which is poor l y understood. Recent re- 



sults ('Oppenheimer et aL 201C 
[Pauchcr-Gigucrc 



. .van de Voort et all I2OIII : 

Kcrcs. fc Mai boilh have demonstrated 



that this can be a significant effect in plausible (though 
somewhat extreme) wind models. Note that although both 
Cwinds and rj arise from winds, the former is a preventive 
feedback parameter, whereas 7] is an ejective feedback pa- 
rameter; the two are not necessarily related in a simple way, 
so we keep them separate. 

Multiplying these terms together, we obtain C(-^haio) 
that is schematically illustrated in Figure [1] It is small at 
low and high halo masses owing to photo-suppression and 
quenching respectively, and approaches unity at intermedi- 
ate halo masses which is where vigorous star formation can 
occur. This depiction is intended only to illustrate broad 
trends, as the actual values and functional forms of the var- 
ious C terms are at best only qualitatively known. The gen- 
eral shape of this curve with a cut off in accretion at low and 
high masses is long known (e.g . Thoul fc Weinberg Il996l : 
iKeres et af] l2005l : ICroton et al.l l2006l ). although the quan- 



titative masses for the cutoffs are debated to this day (for 
recent work on this see iBouche et alll2010l : 1 Cattaneo et al.l 

Iml). 

There may be other sources of energetic preventive feed- 
back that retard accretion such as cosmic rays, stellar winds, 
quasar outflows, local photo-ionisation, and magnetic fields, 
but their importance for global galaxy evolution has not yet 
been firmly established. There may also be subtle "amplifi- 
cation" effects by which two or more preventive mechanisms 
serve to strengthen each other beyond their individual im - 
pact jPawlik fc Schayell2009l : iFinlator. Dave, fc 6zej|201ll ). 
Hence the list of individual C's above is intended to illus- 
trate the sort of physical processes contributing to preven- 
tive feedback, and how they might manifest in the overall 
shape of C- The actual trend of C(Mhaio) may involve more 
complicated and subtle effects than described here. 

Figure [2] illustrates the impact of these various in- 
flow and feedback terms on the speciflc star formation 
rate (sSFR=SFR/M,) as a function of Mhaio, using large- 
scale cosmological h ydrodynamic simulations with various 
outflow models (see IPave. Oppenheimer. fc Finlatoil l201ll . 
for details). The upper left panel shows the case with- 
out outflows. If all gas entering into the halo ended up 
in the ISM, i.e. A/in = Mgrav, then the relation would 
be as shown by the solid line, having a positive slope. 
Even without any feedback, gravitational heating results 
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Figure 2. Specific star formation rate (sSFR=SFR/A/» ) versus 
halo mass M^^^l^ at z = in 48 Mpc/h cosmological hydrody- 
namic simulations (see lDave. Oppenheimer. &: Finlatoi| [2011. for 
details). Each point represents a galaxy, and at a given mass 
these are subdivided into high /gas (blue; > 0.5cr above the me- 
dian), low /gas (red; < — O.Scr below the median), and inter- 
mediate (green). Upper left panel shows the case with no out- 
flows (nw|^_So]M_b]ackJliie_s^ scaling of Mgrav 
from iFakhouri. Ma. fc Bovlan-KilchinI l|2010l '). while dashed line 
shows the effect of gravitational pre ventive feedback (eq. [B] and 
iFaucher-Giguere. Keres. fc Ma|[201ll '). The cyan line shows a run- 
ning median; this curve is reproduced in the other panels for com- 
parison. The upper and lower right panels show models where we 
assume r) = 2 and i^wind = 680 km/s and 340 km/s (cw and sw), 
respectively. T he lower left panel shows momentum-d riven wind 
scalings (vzw; iMurrav. Quataert. fc Thompson! |2005^ ■ in which 
r] oc Afj^^^J^"^ and the outflow speed is oc M^(^^^ (roughly). De- 
viations from the no-wind case represent the influence of other 
terms on Mi^, particularly Mrecyc at high masses and Cwinds 
low masses. 



in a negative slope (daslied line). The results from our 
simulations are in good agreement with Equation [6] (from 
iFaucher-Giguere. Keres. fc MallioT Vi), since the simulations 
themselves are quite similar; the main difference is that ours 
include metal-line cooling, but this is not very i mportant un- 
til h ot gaseous halos form at A/haio ^ 10^^ M© (| Keres et al.l 
12005 : Gabor ct al. 2011) since cold accretion is generally lim- 
ited by the infall time rather than the cooling time. 

The simulations in the two right panels assume 77 = 2. If 
Min is unchanged, then one expects both the SFR and M, to 
be lowered by a factor of three (eq.[5]), resulting in no change 
in sSFR from the no-wind case (the no-wind cyan curve is 
reproduced in all panels for comparison). But clearly there 
is a change, which reflects the impact of outflows on Min. At 
large masses, wind recycling (Mrecyc) returns material to the 
galaxy rapidly once the wind speed (680 km/s in the upper 



right, 34 km/s in the lower right) drops below the escape 
velocity l|Oppenheimer et al. 2010l). Hen ce at large masses 
rj is effectively (|Finlator fc Pavj 120081 ) . and sSFR jumps 
owing to Mrecyc. Fot the case of the slower 340 km/s winds 
(lower right), the jump occurs at a factor of eight lower in 
mass as expected from the factor of two difference in wind 
speeds. Winds also affect sSFR at low masses, where it is 
suppressed relative to no winds. This reflects ("winds, which is 
as expected stronger in the case of the higher wind speed. In 
the lower left we show simulations using momentum-driven 
wind scalings of (approximately) rj oc 

^haio^- This flattens 
sSFR(Miiaio), in addition to exhibiting different behaviors 
for wind recycling and suppression. 

These examples illustrate how star formation rates at a 
given mass in simulations are impacted by the various ejec- 
tive and preventive feedback processes described above. For 
instance, the fact that observations find no sudden increase 
in sSFR at any c haracteristic mass (e.g. ISalim et"aL I l2007l : 
iDaddi et al.ll2007l ) suggests that galaxies do not eject mate- 
rial at a character ist ic win d v elocity. 

iB ouchc ct al' ("20101), iDutton et al] (|2010l ). and 
iKrum holz & Dckcl (20lj) also present empirical models 
based on accretion-driven star formation, fo cusing on the 
form a nd evolution of the sSFR. Interestingly, iBouche et al.l 
(|2010l ) demonstrates that a model in which ( = 1 for 
10" < Mhaio < IO^^-^Mq and zero elsewhere nicely 
reproduces some key observed properties of high-redshift 
star-forming galaxies, including the observed lack of sSFR 
amplitude evolution from 2; ~ 7 — > 2. Simulations, in con- 
trast, predict that sSFR evolution tracks the a ccretion rate , 
and thus continues to rise to high redshifts (Dave 2008); 
analytic models that d o not include such a low-mass cutoff 
show similar behavior jPutton et al.ll2O10l ). If this low-mass 
cutoff for accretion were true, it would suggest that there 
are additional feedback processes affecting small systems 
beyond what is shown in Figure [T] namely, metagalactic 
photo-ionisation. An alternative explanation that does not 
employ a sharp mass cutoff but still reproduces the observed 
sSFR behavior is to invoke an observationally-motivated 
me tallicity-dependent star formation law, as explored 
by iKrumholz fc Dekell HoTl). Pushing such observations 
of sSFR out to higher redshifts and lower masses (e.g. 
with t he Cosmic Assemb l y Near-Infrared Deep Legacy 
Survey; iGrogin et all I2OIII ; iKoekemoer et all I2OIII ) should 
provide interesting constraints on the physical mechanisms 
regulating early, low-mass galaxy growth. 



3 GAS FRACTIONS 



A galaxy's gas fraction is defined here as 



A'4 



1 



Mgas + M. 1 + (tdepsSFR) - 1 



(7) 



where in the second equality we have employed the deple- 
tion time tdcp = Mgas /SFR. We argue below that the latter 
formulation offers the intuitive advantage that it splits /gas 
into a term that is fairly insentive to feedback (tdcp) and a 
term that depends strongly on feedback (sSFR). 

The depletion time measures the timescale over which 
gas, when present in the ISM, gets converted into stars. This 
is expected to be primarily determined by the star formation 
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law, such as the observed lKennicut j l| 19981 ) relation between 
gas surface density ( Sgas) to SFR surface density (Ssfr). In- 
deed, simulations bv lPave. Finlator. Sz Oppenheirneil l|201ll . 
hereafter DFOll; see their Figure 4) assuming a Kennicutt- 
Schmidt law show that tdcp is essentially independent of 
outflows, and scales as 



idop oc tnM' 



(8) 



where tH is the Hubble time. 

We can derive this scaling of tdcp directly from the star 
formation law. The temporal scaling can be most easily un- 
derstood using a formulation of the star formation law given 
by SFR« 0.02Afgas/fdyn, where tdvn is the dynamical time of 
the s t ar formation r egion (e.g. ISili3ll997l : iKrumholz fc TanI 
I2OO7I : iGenzel et all I2OI0I '). This then gives tdep oc tdyn, 
whi ch in a canonic al disk model scales as the Hubble time 
tH (|Mo. Mao, fc White 19981. Meanwhile, the stellar mass 
dependence arises from the Kennicutt law plus the ISM gas 
profile. The Kennicutt law states that Esfr oc (where 
A'' ~ 1.4), from which it is straightforward to show that 
tdep oc T^l-f. In simulations, Egas oc M^^'^{l + zf (DFOll), 
which gives rise to the weak anti-correlation with M* quoted 
above. We caution that these dependences may be somewhat 
different in the real Universe since these simulations lack the 
resolution to properly model the internal structure of galax- 
ies. 

The evolution of /gas depends on the evolution of tdcp 
and sSFR. The former evolves with Ih (e.g. as (1 + 2:)"^'^ in 
the matter-dominated regime), while the latter is generally 
driven by cosmic inflow (eq. [3| which scales as (1 -I- z)'^'^^ 
if driven by gravitional infall. Combining these, galaxy gas 
fractions are predicted to evolve with time roughly as ~ 
(1 -I- z)'^''^^tH (if /gas is not near unity), which increases 
slowly with redshift , qual i tatively consistent with observa- 
tions l|Tacconi et al.ll2010l : ICeach et al.ll201ll '). Hence in the 
equilibrium scenario, galaxy gas fractions represent a compe- 
tition between supply and consumption, such that galaxies 
become less gas-rich with time because the gas supply rate 
OS faster than the gas consumption rate. 



4 METALLICITIES 

The global metallicity within the ISM is given by the enrich- 
ment rate, which is the yield y times SFR, divided by the 
mass inflow rate Mjn t hat must be enriched. As derived in 
iFinlator fc Davj l|2008l ) , if the inflow is pre-enriched there is 
an additional term that depends on az = Zin/ZisM, where 
Zin and ZisM are the metallicities of the inflowing and am- 
bient ISM gas, respectively: 



■^iSM = y 



SFR 

Min 



y 



1 + 77 1 — az 



(9) 



Hence the mass-metallicity relation and its evolution 
are established by a the mass and redshift depen- 
dence of ri and az- In our currently fav ored outflow 
model (|Dave. Qppenheimer. fc Finlatoil(201ll ). rj has a sig- 
nificant mass dependence but little or no redshift depen- 
dence, while az is generally small but has a significant red- 
shift dependence. Equation [5] would then suggest that the 
shape of the mass-metallicity relation is primarily estab- 
lished by ri{M,), while its evolution is driven by az; this was 



demonstrated for simulations in DFOll. Hence in this sce- 
nario, the shape of the mass-metallicity relation is modulated 
by the fraction of inflow that forms stars, whtle its evolution 
IS governed by the enrichment level of inf ailing gas. 

Conspicuously absent in this scenario is any explicit ref- 
erence to potential wells of galaxies, or any consideration of 
outfiow velocities versus escape velocities. These processes 
are canonic ally believed to govern the mass-me tallicity re- 
lation (e.g. iDekel fc Sil"^ Il986l : iTremonti |2004): the 
phrase "metals can more easily escape from the shallower 
potential wells of small galaxies" is oft-repeated. However, 
in our scenario, it is instead the net mass outflow rate that is 
the key determinant of the mass-metallicity relation, and the 
potential well depth is at most only indirectly implicated. 

Since the vast majority o f metals in the IGM are de- 
jpositcd t here by outflows (e.g. lOppenheimer fc Davellioosi : 
[Oppcnhci mer et al.ll201lh . the infalling gas metallicity is a 
direct measure of Mrecyc- Zi-n is given by the metal mass 
arriving in the form of recycled winds, divided by the total 
mass inflow rate, i.e. 



Zin — Zy^ 



Mr, 



(10) 



where the denominator is Min from Equation [5] 

Under the typical case of highly mass- loaded outflows, 
the outfiowing metallicity must be similar to the ambi- 
ent ISM metallicity, i.e. Zout ~ ^ism. Furthermore, since 
galaxies evolve slowly in metallicity (e.g. iBrooks et al ] |2007l . 
DFO ll) and wind r ecycling times are typically of order a 
Gyr (jOppenheimer e t al. 2010), the galaxy metallicity has 
probably not evolved strongly from when the gas was ejected 
to when it is being re-accreted, and hence Zrccyc ~ ZoMt. 
Substituting Zrccyc ~ ■Z'isM into Equation [10] and solving for 



Mr, 



yields 



az 

Mrocyc = CMgrav. (H) 

1 — az 

This relates the mass recycling term in the inflow equation 
(eq. [2| to the metallicity infalling into the ISM. The advan- 
tage of formulating recycling in terms of qz is that it is in 
principle an observable quantity via absorption or emission 
measures in the outskirts of galaxies. In contrast, Mrecyc is 
not directly measurable since it is not clear how to distin- 
guish recycled wind inflow from other inflow, or even how 
to measure galaxy inflow rates at all. Note that since galaxy 
metallicities evolve slowly upwards with time. Equation 1111 
will tend to slightly underestimate Mrecyc for a given az. 
The ejection of winds from one galaxy (typically a satellite) 
being accreted onto another (typically the associated cen- 
tral) could also affect az, which would also cause an under- 
estimate in Mrecyc slucc the satellites are generally smaller 
and hence lower metallicity. 

5 THE EQUILIBRIUM RELATIONS & 
IMPLICATIONS 

We can substitute Equation [TT] into Equation [S] to obtain 

(l + 7))(l-az) 

This is the key equation that delineates how galaxy star 
formation rates are governed by accretion and feedback 
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processes, i.e. baryon cycling. This equation, together with 
Equations [7] and [9l represent the equilibrium relations that 
govern the steUar, gas, and metal content of galaxies across 
cosmic time. Galaxies will tend to lie around these relations 
owing to a balance of inflow, outflow, and star formation. 

The equilibrium relations depend on three parameters: 
rj, and az, representing ejective feedback (i.e. outflows), 
preventive feedback, and wind recycling. Additionally, the 
star formation law governs tdop, Afgrav is set by cosmology, 
and y is set by nucleosynthetic processes. Assuming those 
are well-established, the mass and redshift (and possibly en- 
vironmental) dependence of rj, and az govern the evolu- 
tion of the global SFR, /gas, and Zism of galaxies. Note that 
since the mass and redshift dependence of these parameters 
are not fully known, the actual number of free parameters 
can be significantly larger than three. 

There are many possible ways to chara cterize simula- 
tion results into an analytic formalism (e.g. iNeistein et al.l 
One virtue of our particular parameterization is that 
the parameters involved are, at least in principle, directly ob- 
servable. This provides an optimally direct connection from 
observations to constraints on galaxy formation models. Un- 
fortunately, measuring these parameters is challenging, but 
preliminary constraints have already been obtained. 

For instance, rj has been constrained in hi gh-g galaxies 
to have a value of orde r unity or more (e.g. ISteidel et al.l 
I2OIOI : iGenzeT et al.l[201ll ). az can be constrained by exam- 
ining metallicities in the outskirts of low-2 galaxies (e.g. 
iBresolin et~al]|2009l : [Moran et al.ll201ll '). Constraining C by 
direct observations would require an accurate census of all 
halo gas which is highly challenging, but aside from ("winds, 
its main terms can be constrained using a combination of 
relatively straightforward numerical work ((photo and Cgrav) 
and empirical arguments ((quench). 

The equilibrium relations have some interesting impli- 
cations for the behavior of SFR, /gas, and Z. For instance, 
hydrodynamic simulations indicate that the star formation 
history of galaxies is insen s itive to the assume d star forma- 
tion law ijKatz et al.lll996l : ISchave et al.ll2010h . This seems 
paradoxical at first, but is straightforwardly seen from Equa- 
tion [T2j since there is no dependence here (or in the metal- 
licity equation) on the star formation law. The star forma- 
tion law only affects the gas fractions, via f dep. This can 
be re garded as a self-regulation mechanism (|Schave et al.l 
in which gas collects in galaxies as required in or- 
der to achieve the star formation rate set by the balance of 
inflows and outflows. 

Another straightfoward prediction of the equilibrium 
relations is that if one desires the mass-metallicity re- 
lation to scale as Z (X mI ^ ^ at sma l l ma sses as 
observed (|Tremonti et al.l |2004| : iLee et all 120061 ') . then 
Equation [9] directly implies rj oc M^^^^ (assum- 
ing az <C 1), roughly as expected for momen tum- 
driven winds (|Murrav. Quataert. fc ThompsonI |2005| ). In- 
deed, simulations assuming such a scaling appear to pro- 
vide a good match to mas s-metallicity relation observa- 
tions l|Finlator fc Davell2008l . DFOll). 

It is instructive to combine Equations [9] and [12] to give 



SFR 



i.e., the fraction of gravitational infall into a halo that ends 
up forming into starfl, while the second ratio quantifies the 
metal retention fraction within galaxies. If the halo mass 
and metal yield can be determined, measuring SFR/Zism 
provides a quantitative constraint on (. This can be done at 
least at 2 ~ with existing data from e.g. the Sloan Digital 
Sky Survey. 



6 A SAMPLE EQUILIBRIUM MODEL 

The equilibrium model can be used to quickly explore pa- 
rameter space and obtain intuition about the governing 
physics for galaxy properties of interest. We illustrate this 
here by presenting results for the evolution of galaxies in a 
full equilibrium model. 

Figure |3] shows the evolution of the SFR, halo SFE, Z, 
and /gas for four galaxies spanning the indicated range of 
final [z = 0) halo masses. These are computed using Equa- 
tions [71 [9l and 1121 tracking the stellar and halo mass growth 
starting at an early epoch when the halo is at the photo- 
suppression mass. We take Mgrav from Equation[3l parame- 

/q)""'^ as discussed in [jS] and 



Mgrav ZlSlA 



(13) 



The first ratio is the halo star formation efficiency (SFE), 



grav 

terize tdep = Q At h [M, /IQ^° Mq)-°-^ 
take az = (0.5 - 0.lz){M, Mq)°-'^^ (with az ^ 0) as 
a crude parameterization of simulation results from DFOll. 
We choose ( as described in Figure [1] and define rj as indi- 
cated in the upper right of the Figure: solid lines approx- 
imately represent momentum-driven wind scalings, while 
dotted lines represent energy-conserving wind scalings. We 
also include instantaneous recy cling of 18% of s tar formation 
back into gas as expected for a lChabrieil (|2003l ) IMF, but do 
not include further stellar mass loss. We reiterate that these 
"base model" parameter choices are at some level arbitrary, 
and are intended only to illustrate how parameter variations 
influence observables. 

The green line represents a Milky Way-sized halo of 
IO^^Mq. At z = 0, it has SFRfa 2.5MQ/yr, 5 x 

lO^^M©, /gas ~ 0.1, and Z ~ Zq, in fair agreement with 
measured values and showing that our parameter choices 
are reasonable. Larger galaxies form stars more vigorously 
at earlier epochs and for shorter intervals, which is quali- 
tatively similar to the beh avior in the empiric al "staged" 
galaxy formation model of iNoeske et al.l |2003). The peak 
SFRs are 50-100 Mq/jt at z ~ 2 — 3, which is lower than 
the observed values for the largest main sequence galaxies at 
that epoch by a factor of a f ew, reiterating the issue noted 
in simulations bv lDavj (|2008l ) that observed galaxy SFRs at 
that epoch approach or exceed their cosmic accretion rate; 
the resolution to this quandary remains un clear (see e.g. 
iBouche et al]|2010l : IXrumholz fc Dekelll201ll . for ideas). 

The halo SFE is plotted in the second panel. For star- 
forming galaxies, this efficiency is roughly one-third to one- 
half over most of cosmic time. Higher mass halos show a 
marked drop in efficiency once they grow above the quench- 
ing mass, here assumed to be 2 x 10^^ Mq at all epochs. 
Going to a steeper scaling of r)(Mh) oc M^"^^^ (dotted 



^ This is distinguished from the ISM SFE, which is how much gas 
is converted into stars over some characteristic galaxy timescale, 
or the "cosmological" SFE, which is the galaxy stellar mass 
divided by the cosmologically-expected halo baryon mass (i.e. 

Af,//i,Afhalo)- 
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Figure 3. Equilibrium model evolution of four galaxies that 
have final z = halo masses of IO^Mq, 10^^ Mq, IO^^Mq, 
and 10^* Mq (blue, green, magenta, and red lines, respectively). 
Panels from top to bottom show the star formation rate, cos- 
mic star formation efficiency defined as the fraction of baryons 
entering the halo that form into stars, gas fractions, and metal- 
licities. Solid lines show r] = {Mf^/lO^^ Mq)~^^^ , and dotted lines 
show 77 = (Mfi/1O^^M0)~^/^. The star formation histories peak 
earlier and at higher SFRs for larger galaxies, r eminiscent of 
the em pirical "staged galaxy formation" model of iNoeske et alj 
Star formation efficiencies are typically ~ 1/2 — 1/3 for 
star-forming galaxies, and drop quickly about the quenching halo 
mass of 10^^ '^Mq for more massive systems. Galaxies self-enrich 
to above one-tenth solar very early on, and then evolve slowly 
in m etallicity. Obse r ved g as fractions for fairly mas sive galaxies 
from Tacconi et al. I ll2010f) at 2 = 1.2 and 2.2 and iGeach et al.l 



(|201l|) at z = 0.4 are indicated, to be compared with predictions 
for ~ 10^^~ ^^Mq halos; the model predictions are generally too 
low. The dashed lines show the experiment of adding a metallicity 
dependence of oc to t^cp when sub-solar, which keeps early 

galaxies gas-rich and improves agreement with data. 



evolve much (DFOll). A steeper scaling of r){Mh) strongly 
steepens the mass-metallicity relation in smaller galaxies at 
all redshifts, as expected from Equation!^ 

/gas drops slowly with time since the accretion rate 
drops faster than the consumption rate, as discussed in ^ 
Here we include some observation al comparisons (large sym- 
bols), from CO me asur ements by|Tacconi et al.l l|2010l ) and 
iGeach et all l|201lh . As lXacconi et all l|201Cl ) noted, the rel- 
atively slow evolution requires continual replenishment, as 
our model naturally predicts. However, the z ~ 1—2 gas frac- 
tions are higher than the model predictions, because the ob- 
served galaxies are best compared to the green and magenta 
lines representing final halo masses of ^ 10^^ — IO'^^Mq. 
The equilibrium model can quickly test whether parameter 
variations can reconcile this discrepancy. 

According to Equation [T] higher /gas values can be 
achieved by raising sSFR and/or tdcp. Raising sSFR e.g. 
by lowering r] would concurrently alleviate the discrepancy 
with observed z ^ 2 SFRs, but would likely overproduce 
stars globally. Delaying star formation by accumulating gas 
could explain the discrepancy, though it seems somewhat 
contrived to have galaxies at all masses suddenly start vig- 
orously consuming gas around 2 ~ 2. A steeper scaling of 
'r]{Mh) works in this direction, but the dotted lines show 
that this only helps marginally; a redshift-dependent r] that 
is higher at high-z may work better, but the physical mo- 
tivation is not obvious. Another avenue is to raise tdcp by 
appealing to the obs erved lower gas consu mption rate in low- 
metallicity systems (|Bolatto et al.|[201ll l. An illustration of 
this is shown as the dashed lines in Figure [3l where we add 
a dependence of to tdop when below solar metallicity. 
Within this model, we find that this steep metallicity depen- 
dence is required in order to sufficiently raise gas fractions at 
high-z; however, it is not physically-motivated, and does not 
necessarily reflect recent ly proposed metallicity-d ependent 
star formation laws (e.g. iKrumholz fc Dekellboill) . The de- 
lay in star formation results in higher gas fractions at 2 ~ 2 
in better agreement with data, but at lower redshifts when 
metallicities approach so lar this makes less difference. An- 
other avenue explored bv ,Bouche e t al. (2010) is to postulate 
that star formation can only occur in halos above 10^^ Mq, 
which has the advantage of producing higher has fractions 
in smaller systems at earlier epochs; they demonstrate that 
this can broadly match the observed evolution of galaxy gas 
fractions. 

This exercise illustrates how an equilibrium model can 
be utilized to gain intuition about various physical processes 
in galaxy evolution. The relatively small number of param- 
eters and their direct connection to observable physical pro- 
cesses makes this type of model easier to interpret than mod- 
ern simulations or semi-analytic models. 



lines) more strongly suppresses low-mass galaxy growth, 
which is favored in semi-analytic models to reproduce the 
observed faint-end slo pe of the stellar mass function (e.g. 
ISomerville eral]|2008l ). 

Metallicity evolution is very rapid in all galaxies early 
on, and galaxies typically enrich to > O.lZp) within 
the flrst Gyr, as seen in simulation s jPave et al.l I2OO6I : 
iFinlator. Oppenheimer. fc Dav4 I2OIII ). Evolution is slow 
thereafter. Galaxies tend to follow parallel tracks, show- 
ing that the shape of the mass-metallicity relation does not 



7 DEPARTURES FROM EQUILIBRIUM 

Stochastic variations in A/in, including mergers, can cause 
departures from equilibrium. This generates scatter about 
the equilibrium relations. Generically, departures from equi- 
librium tend to return galaxies towards equilibrium. This 
self-regulating behavior is why we dub this scenario the equi- 
librium model. 

To illustrate this behavior, consider a galaxy experi- 
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encing an upward fluctation in Min. Its /gas increases, and 
perhaps M« as well if there is a small galaxy accompanying 
the infall, while Zism decreases since either fresh infall or 
a lower-mass galaxy will have lower metallicity. This moves 
the galaxy off the equilibrium relations. The increased gas 
content immediately stimulates more vigorous star forma- 
tion, which over time enriches the galaxy as it consumes the 
excess gas. This then lowers its gas content and increases its 
metallicity, moving the galaxy back towards equilibrium. 

Conversely, a temporary lull in accretion will make a 
galaxy more gas-poor and metal-rich as it consumes its 
existing gas. Quantitatively, its metallicity will evolve fol- 
lowing the relation A Zism ~ yAAf,/Afgas (see eq. 19 of 
iFinlator fc Davell2008l . in the case of zero accretion), which 
is steeper than the MZR and hence will move the galaxy 
above the MZR. In time, hierarchical growth will bring in 
fresh gas, lowering the metallicity and increasing the gas 
content to return the galaxy towards equilibrium. In this 
way, galaxies oscillate around the equilibrium relations, con- 
stantly being perturbed from them and driven back owing 
to fluctuations in infall. 

An inevitable prediction of this scenario is that 
departures from the equilibrium relations will correlate 
with star formation rate, gas fraction, and metallicity. 
From the above scenarios, one can see that at a given 
mass, galaxies that are gas-rich (gas-poor) and metal- 
poor (metal-rich) will have higher (lower) star forma- 
tion rates. T hese trends are qua li tatively consisten t with 
observations ([Ellison et al.l 1200^: iLara-Lopez et al.l l2010l : 



iMannucci et al.ll2010l : IPeeples fc Shankaj|201of ) 



This trend is illustrated for gas fractions and star for- 
mation rates by the coloured points in Figure [51 showing 
that at a given mass, high SFR and high /gas go hand in 
hand. Figure 1 of DFOll analogously shows that high SFR 
accompanies low Zgas at a given M, . Note that these second- 
parameter trends do not arise from outflows, being present 
even in simulations without winds. Instead, it is a direct and 
unavoidable consequence of equilibrium, and results from 
galaxies' self-regulating response to fluctuations in inflow 
rather than any feedback process. 

Quantitatively, the scatter around the equilibrium rela- 
tions depends on how quickly galaxies can return to equi- 
librium after being perturbed. To return, there must be 
sufficient infall to re-equilibrate the galaxy. The timescale 
for t his to happen can b e quantified by the dilution 
time (jFinlator fc DavjbOOsI '): 



idil 



Mg, 



(1 + 77) 



.lAfg 



M, 



(l + 77)~'tdcp. (14) 



If the dilution time is small compared to the in- 
flow fl uctuation timescale, t hen the scatter will be 
smaU (|Finlator fc Davd 120081 '). The inflow fluctuation 
timescale likely depends on mass and environment, with 
small galaxies generally suffering (relatively) larger pertur- 
bations. Equation 1141 shows that the dilution time depends 
on rj and tdcp, and hence observations of the scatter ver- 
sus M« provides an independent constraint on r]{M-t) and 

WLM*1 

lOutton et all (|201(]| ) argues, based on an analytic model 
of accretion-driven galaxy formation similar to this one, 
that the scatter in the observed Af, — SFR relation must 
be driven by fluctuations in inflow, since assuming zero 



scatter in the relation between Afin and Afhaio results in 
a relation that is too tight compared to observations. How- 
ever, hydrodynamic simulations that implictly include in- 
flow fluctations also yield a small sca tter (e.g. JDavc 200^; 
IFinlator. Oppenheimer. fc Pavel 1201 ll *) . Therefore it is not 
solely inflow fluctuations that govern the scatter, it is the 
more complex relationship between inflow fluctuations and 
the dilution time. Equation [14] suggests that this, in turn, 
depends on the outflow rate {rf) and gas consumption rate 
(tdcp). The fact that feedback regularizes galaxy properties 
may have other effects on galaxy properties. For instance, it 
has long been suggested that the low scatter in the TuUy- 
Fisher relation arises owing to feedback processes, since fluc- 
tuations in halo growth alone w ould naively predict a scat ter 
that is large compared to data (|Eisenstein fc Loeblll996l ). 

Satellite galaxies lie permanently off the equilibrium re- 
lations, because the inflowing filaments bypass them and 
flow to the centers of halos. Hence they are expected 
to end up with lower gas content and higher metal- 
licities than centrals of the same m ass; this is as ob- 
served l|Peeples. Pogge, fc Stane^ 120091 '). How far they lie 
off the equilibrium relations depends on their gas reservoir 
at the time they are cut off from their supply, which sim- 
ulations indicate is typ ically ^ 1 Gyr af ter falling into a 
hot gas-dominated halo (|Simhaet al.li2009l ). Note that some 
satellite galaxies fall in along the inflowing filaments, and in 
those cases they are not bypassed since they are actually 
part of the inflow. However, these satellites are expected to 
quickly merge into the central galaxy, and so will not typi- 
cally end up as part of the long-lived satellite population. 

Major mergers are another population lying far out 
of equilibrium. It is not merely that such systems rep- 
resent a particularly large inflow perturbation, it is that 
the induce d torques drive gas fiows that fuel central star 
formation (|Mihos fc Hernguistl 1199^ ) making cosmological 
inflow mostly irrelevant during the merger event. In a 
global context, major mergers are seen to be responsible for 
only a small frac tion of overall cosmic star formation (e.g. 
^Jog ee et al.ll200^ '). although they may be where muc h of the 
central black hole mass growth occurs (e.g. Di Matte o~et al.l 
l2005h . Overall, the equilibrium model is broadly valid for 
central galaxies th at are quiescently for ming stars, i.e. main- 
sequence galaxies (|Noeske et al.|[2007t ). that dominate cos- 
mic star formation. 



8 BEFORE EQUILIBRIUM: THE GAS 
ACCUMULATION PHASE 

Since the infall rate has a steeper redshift dependence 
than the consumption rate (e.g. ~ (1 -f z)^'^^ vs. t^), 
at sufficiently early epochs galaxies will be in a gas accu- 
mulation phase during which galaxies cannot process gas 
into stars as fast as th ey receive it (|Bouche et all l2010l : 
iKrumholz fc Dekej|201ll l . Only when consumption can keep 
up with supply will equilibrium be achieved. During the gas 
accumulation phase, gas fractions are expected to be higher 
and metallicities lower than predicted by equilibrium. 

'We can estimate the redshift Zoq where gas accumula- 
tion ends and equilibrium is attained. Star formation must 
be able to occur fast enough to satisfy SFR= A'/in/(l -f rf) 
(eq. [2]). We assume that at early times, ~ 1 and az ~ 0, 
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Figure 4. Evolution of the ratio of the infall term (eq. [2]| versus 
the SFR term (SFR=0.02Afgas/tdyn), for different choices for r] 
and halo mass (solid and dashed lines represent Af;, = 10^^ and 
VS^^Mq, respectively). This calculation assumes /gas = 0.5 for il- 
lustration purposes, which is typical at high- 2: but in detail varies 
with r\ and My^. With no outflows (r; = 0; blue lines), equilibrium 
for massive galaxies occurs at 2 ~ 2, while less massive galax- 
ies equilibrate earlier. With modest outflows (r; = 1; red lines), 
equilibrium is already pushed back into the early Universe. The 
green lines assume momentum-driven wind scalings for T) from 
lOppenheimer fc Dav e (2008), which pushes equlibrium to 2 ~ 5 
for massive galaxies, and much earlier for smaller ones. 



so that Min ~ Mgrav (eq. [3}. We further assume that 

10*(1 -I- zY^I"^ yr is the 
- /gas/fcMhaio, we obtain 



SFR=0.02Mgas/tdyii, where tdyn - 
disk dynamical time. Taking Mgas 
an equation for Zoq: 



1 +2o 



[^/-<-')1"'(t^) 



(15) 



We illustrate these trends in F igure [H For sizeable high- 
z galaxies with /gas ~ 0.4 — 0.5 (jTacconi et al.l[201ol ) and 
no outflow s (77 = 0), we obtain Zeg ~ 1.5 — 2.5, in agree- 
ment with iKrumholz fc Dekell | 201 ll ). But the superlinear 
dependence on /gas and (1 -I- r]) means that the results are 
quite sensitive to these values. For instance, even a modest 
outflow rate of = 1 yields ~ 5 — 7. Effectively, out- 
flows lower the amount of inflow that needs to be processed 
into stars, allowing for an earlier equilibration epoch. There 
is a weak halo mass dependence as well such that smaller 
galaxies equilibrate earlier. Note that this derivation of Zoq 
depends on the star formation law: if the star formation 
law was different in early galaxies owing e.g. to metallicity 
effects, then this would also impact when galaxies achieve 
equilibrium. As such, a precise prediction of Zeq is sensitive 
to poorly known factors. Nevertheless, with realistic out- 
flows, it is likely that Zeq S> 2, and hence galaxies live in 
equilibrium over the vast majority of cosmic time. 

Observationally, IPapovich et al.l (|2010l ) used the star 
formation rates, masses, and (estimated) gas contents of 
high-z Lyman break galaxies to infer that gas accumulation 
occurs down to Zeq ~ 4, after which accretion and star for- 
mation track each other as expected in equilibrium. Hence 
there is some direct empirical support for an early gas ac- 
cumulation epoch. Constraining jZcq more precisely will pro- 



vide quantitative constraints on gas processing rates in early 
galaxies. 



9 SUMMARY AND DISCUSSION 

We have presented a simple formalism for understanding 
the evolution of the stellar, gaseous, and metal content of 
galaxies, inspired by intuition gained from cosmological hy- 
drodynamic simulations. This formalism is encapsulated by 
the equilibrium relations: 



SFR = 



.^ISM ~ 



(l + 77)(l-az)' 
1 

1 + (tdopsSFR)-i 

y ^_ 



(16) 
(17) 
(18) 



1 + 77 1 

These relations are established by a balance between inflows 
and outflows, and evolve on cosmological timescales over 
which inflow and outflow rates slowly vary. They are pri- 
marily governed by three baryon cycling parameters that de- 
scribe ejective feedback (r\) , preventive feedback (C^ , and the 
re-accretion of ejected material (az); each of these param- 
eters is in principle directly observable, but at present has 
poorly known dependences on mass and redshift (and per- 
haps other properties). Additionally, they depend on Mgrav, 
the gravitational infall rate of baryons into the halo set by 
ACDM, tdcp, the ISM gas depletion time, and y, the metal 
yield. These relations capture, to flrst order, the behavior 
of galaxies in modern hydrodynamic simulations that in- 
corporate these processes dynamically within a hierarchical 
structure formation scenario. 

The equilibrium model is broadly valid for quiescently 
star-forming central galaxies, at epochs where star formation 
is able to keep up with inflow (i.e. past the gas accumula- 
tion epoch; Figure [J), and when averaged over timescales 
longer than stochastic fluctuations in the inflow rate. On 
shorter timescales, galaxies oscillate around the equilibrium 
relations such that more (less) rapidly star-forming galax- 
ies at a given mass having higher (lower) gas fractions and 
lower (higher) metallicities. The scatter about the relations 
is governed by a competition between the dilution time 
tdii = (1 -|- ri)~^tdep and the inflow stochasticity timescale. 
This model is not valid for satellite galaxies disconnected 
from feeding fllaments, or for galaxies undergoing a major 
merger where gas feeding is temporarily driven by internal 
dynamical processes. Nonetheless, observations indicate that 
quiescently star-forming galaxies along the so-called galaxy 
main sequence domin ate cosmic star formation at all epochs 
where measured (e.g. iNoeske et al.l l200"7l : iRodhigiero et al.l 
120111 ). and hence this model describes how the bulk (but 
not all) of the stars in the Universe formed. 

Most current galaxy formation models, both hydrody- 
namic and semi-analytic, already include inflow and outflow 
processes within growing large-scale structure. Hence there 
is no new physics in the equilibrium model. What is no- 
table is not what this scenario contains, but rather what it 
doesn't contain. In particular, there is no explicit mention of 
mergers, disks, environment, cooling radii, or virial radii — 
all central elements in the canonical scenario for galaxy for- 
mation. Such elements are automatically accounted for in 
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hydrodynamic simulations, which form disks, merge them, 
and imphctly include environmental effects within growing 
large-scale structure. Yet the equilibrium relations well de- 
scribe such simulations without reference to these elements, 
suggesting that they are not of primary importance for the 
evolution of galaxies' SFR, /gas, and Z. 

In a broader context, the usefulness of the equilibrium 
scenario is that it lays bare the overwhelming complexity 
of modern galaxy formation models, and isolates those as- 
pects that are critical for governing global galaxy evolution, 
thereby providing a simpler intuitive view for how galaxies 
grow. Although the canonical "halo-merger" view of disks 
cooling within halos and merging to drive galaxy evolution 
is not incorrect (i.e. these processes do happen), such a view 
obfuscates the primary driver of global galaxy evolution, 
namely the balance between inflows, outflows, and star for- 
mation. Indeed, the very notion of a galaxy halo, which is 
central to the classical view of galaxy formation, is only a 
second-order effect in the equilibrium scenario: The equilib- 
rium relations are driven by the total inflow rate, and the 
"lumpiness" of that inflow owing to individual halos merely 
manifests as scatter around these relations. 

We have argued that the equilibrium model provides a 
reasonable description of sophisticated simulations, but this 
in no way guarantees that it accurately describes the real 
Universe. It is encouraging that certain unavoidable predic- 
tions such as continual gas replenishment and the second- 
parameter dependence of the mass-metallicity relation seem 
to be in broad agreement with observations. But much work 
remains to be done in order to fully test this scenario. In 
particular, the equilibrium model centrally invokes a con- 
tinual cycle of baryons flowing in and out of galaxies as 
a key moderator of galaxy evolution. But direct observa- 
tional evi dence for such processes is currently scant (see e.g. 
iRubin et al. 2011). Critically testing and constraining these 
baryon cycling processes, particularly within circum-galactic 
gas where such processes are likely to be most prominent, 
will be a key contribution from upcoming multi-wavelength 
observational facilities. 

The equilibrium model provides a re-parameterized 
framework for understanding certain key governing aspects 
of galaxy evolution, but is far from a full solution to the 
problem. To fully solve galaxy evolution, we must at min- 
imum understand the physics that governs r;, i^, and az, 
which will require concerted efforts on both observational 
and theoretical fronts. Furthermore, halo accretion rates, 
metal yields, and the star formation law remain uncertain, 
particularly in regimes such as small low-metallicity galax- 
ies. Perhaps most importantly, this scenario in its current 
form explicitly does not address many interesting aspects of 
galaxy evolution such as the establishment of the Hubble 
sequence and the growth of central black holes. It also does 
not include processes that may be central to the evolution 
of certain classes of galaxies; for instance, it does not ac- 
count for stellar ("dry") mergers which are important for 
the late-time growth of large passive systems. Hence much 
work remains to be done in order to comprehensively under- 
stand how galaxies evolve from primordial fluctuations into 
their present state. It is hoped that the equilibrium model 
provides, in its simplicity, a useful guiding framework for 
understanding the increasingly complex problem of galaxy 
evolution. 
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